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Rotational Joint Limits: 1 DoF 



Image: Autodesk, Creative Commons 
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Rotational Joint Limits: 3 DoFs 
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A 


Parameterize 3-DoF Rotations 

• Ideally, we want a 3D parameter space 
that is free of singularities, ... 

• ... in which range bounds for shoulder and 
hip joints can be intuitively expressed. 

• If only 3D rotations would play so nicely... 
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Euler Angles 

• Parameterize 3D rotation by angles of 
rotation about three predefined axes. 

• Choice of axes is arbitrary, as long as no 
two consecutive axes are the same (for 
example, XYZ, ZYX , XZX, YXY, ...) 

• Limit each angle independently. 
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Euler Angles (Cont'd) 




GDC 


GAME DEVELOPERS CONFERENCE March 14-18, 201B • Expo: March 16-18, 2016 #G0C16 


Euler Angles (Cont'd) 

• Euler angles space has singularities due 
to collapsing axes (aka gimbal lock). 

• Only suitable for joints formed by nested 
gimbals. 

• Totally unsuitable for shoulder and hip 
joints. 
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Euler's Rotation Theorem 


"Any orientation of a 3D object 
can be reached from an initial 
orientation by executing a single 
rotation about a suitable 3D axis. 

(Leonard Euler, 1775) 
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Axis-Angle Parameterization 

• Axis is represented by a normalized 3D 
vector (3 parameters, 2 DoFs!). 

• Zero-angle rotations form a singularity 
(axis is arbitrary). 

• Hard to express joint limits. 
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Exponential Map Parameterization 

• Scale axis by angle to form a 3D vector 
with three independent parameters. 

• Zero-angle rotation is represented 
uniquely by the zero vector. 

• Still has singularities for angles that are 
multiples of 2n (360°). 
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Exponential Map (Cont'd) 

• Limit angle range to [0, 2n). This clears 
out all singularities. 

• We still have a double covering. Rotating 
with angle 0 about axis u results in the 
same orientation as rotating with angle 
2n -6 about axis -u. 
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Exponential Map (Cont'd) 

• Limiting the angle range to [0,7r] restricts 
the double covering to angles of n. 

• Parameterization space is a 3D ball with 
radius n. 

• Admissible orientations form a volume 
inside the 3D ball. 
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Quaternions 

• Quaternions extend complex numbers 


q = w + xi + yj + zk 

where w, x, y and z are real numbers 

• w is the real or scalar part, and 

• (x, y, z) is the imaginary or vector part. 

i it. 
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Quaternions (cont'd) 

• Quaternions behave as 4D vectors w.r.t. 
addition and scaling. 

• In multiplications, the imaginary units 
resolve as: i 2 = j 2 = k 2 = ijk = -1 

• In scalar-vector notation, multiplication is 

given by: \w lr v 2 ] = 

[l W 1 W 2 ~ Vi • V 2 W{V 2 + W 2 V 1 + v x x v 2 ] 
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Unit Quaternions 

• Unit quaternions (points on sphere in 4D) 
form a multiplicative subgroup. 

• A rotation with angle 6 about unit vector 
u is represented by unit quaternion 



Unit Quaternion Parameterization 

• Unit quaternions are parameterized by 
four dependent parameters (3 DoFs). 

• 3D orientations are doubly covered: q 
and -q represent the same orientation. 

• Yet, for procedural animation quaternions 
are generally the best choice. 
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Who Needs the Scalar Part? 

• The scalar part w, less a sign, can be 
found given the vector part v, since 

W = ±Vl — v ■ V 

• Any orientation can be represented by 
unit quaternion with nonnegative w. 
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Quaternion Vectors 

• Quaternion vectors parameterize 
orientations using three independent 
parameters. 

• Zero vector is zero-angle rotation. 

• Parameterization space is a unit ball. 

• Only rotations over n (unit vectors) are 
doubly covered. 
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Quaternion Vectors (cont'd) 

• Quaternion vectors map one-to-one to 
exponential map vectors. 

• For a quaternion vector v, the 
corresponding exponential-map vector is 

2 arcsin(||v||) 

n i v 

v 
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Quaternions to Exponential Map 
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Quaternion Vectors Demo 


& IK Testbed - □ X 

File 

: d 



Ready. 
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Swing-Twist Decomposition 

• Decompose rotation into a swing and 
twist component... 


*1 — Qswing Qtwist 

• ...and, limit each component 
independently. 
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Swing-Twist Decomposition 
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Swing-Twist Decomposition 

• For q = w + xi + yj + zk we find that 

wy — xz wz + xy 

Qswing = s "* “ j ^ 

_ W X . 

Qtwist — 7 + 7 1 ' 


s = Vw 2 + x 2 


where 
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Clamp Swing to Elliptical Disk 


• Quaternion Vector vs. Exponential Map 
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Clamp Swing to Elliptical Disk 

• Map 2D points ( ,k) outside the ellipse to 
their closest point on the ellipse. 

• The outside point lies on the line that is 
normal to the ellipse and passes through 
its closest point. 
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Clamp Swing to Elliptical Disk 
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Clamp Swing to Elliptical Disk 


• A closed-form solution requires solving a 
quartic (4 th order) polynomial. 

• Root finding using Newton-Raphson is 
more accurate, and potentially faster. 

• Incremental error correction is better 
suited for smooth animation. 
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Newton-Raphson 
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Volume Limits 

• Generalization of clamp to ellipsoid or 
elliptic cylinder is straightforward. 

• Clamping to other convex shapes can be 
done using Gilbert-Johnson-Keerthi 
(GJK). 
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Thank You! 

Check me out on 
• Web: www.dtecta.com 


Twitter: @dtecta 

Sample code available in MoTo: 
https://Qithub.com/dtecta/motion-toolkit 


